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Abstract. Consider two insurance companies (or two branches of the 
same company) that have the same claims and they divide premia in 
some specified proportions. We model the occurrence of claims accord- 
ing to a Poisson process. The ruin is achieved if the corresponding two- 
dimensional risk process first leave the positive quadrant. We will consider 
two scenarios of controlled process: refraction and the impulse control. In 
first one the dividends are payed out when two-dimensional risk process 
exit fixed region. In the second scenario whenever process hits horizontal 
line, the risk process is reduced by paying dividend to some fixed point 
in the positive quadrant and waits there for the first claim to arrive. In 
both models we calculate discounted cumulative dividend payments until 
the ruin time. This paper is an attempt at understanding the effect of 
dependencies of two portfolios on the the join optimal strategy of paying 
dividends. For the proportional reinsurance we can observe for example 
the interesting phenomenon that is dependence of choice of the optimal 
barrier on the initial reserves. This is a contrast to the one-dimensional 
Cramer-Lundberg model where the optimal choice of barrier among opti- 
mal barrier strategies is uniform for all initial reserves. 
Keywords: dividend, two-dimensional risk process, proportional reinsur- 
ance 

MSG 2000: 60J99, 93E20, 60G51 

1 Introduction 

In collective risk theory the reserves process X of an insurance company is 
modeled as 

X{t)=u + ct- S{t), (1) 
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where u > denotes the initial reserve, 



Nt 



S{t)^J2^, (2) 



1=1 

is a compound Poisson. We assume that Ui {i = 1,2, ...) are i.i.d. distributed 
claims (with the distribution function F). The arrival process is a homogeneous 
Poisson process Nt with intensity A. The premium income is modeled by a 
constant premium density c and the net profit condition is then c > — , where 

i = nui] < oo. 

Recently, several authors have studied extensions of classical risk theory towards 
a multidimensional reserves model ([1]) where X(t), x, c and S{t) are vectors, with 
possible dependence between the components of S(t). Indeed, the assumption 
of independence of risks may easily fail, for example in the case of reinsurance, 
when incoming claims have an impact on both insuring companies at the same 
time. In general, one can also consider situations where each claim event might 
induce more than one type of claim in an umbrella policy (see Sundt |41|). 
For some recent papers considering dependent risks, see Dhaene and Goovaerts 
[Hill], Goovaerts and Dhaene [21], Miiller [MllSZj, Denuit et al. [TJ, Amba- 
gaspitiya [5], Dhaene and Denuit [T7], Hu and Wu and Chan et al. [T^. 

In this paper we consider a particular two-dimensional risk model in which 
two companies split the amount they pay out of each claim in fixed proportions 
(for simplicity we assume that they are equal), and receive premiums at rates 
ci and C2, respectively. That is. 



Nt Nt \ 

\ 1=1 1=1 / 

Note that ([T]) models proportional reinsurance dependence. The same two- 
dimensional risk process was already considered in Avram et al. [6l [7l [8]. In 
Avram et al. [7j the asymptotics of the ruin probability when the initial reserves 
of both companies tend to infinity under a Cramer light-tail assumption on the 
claim size distribution is analyzed. In Avram et al. [HI 13 foi' the simplest 
particular case of exponential claims the Laplace transforms of few perpetual 
ruin probabilities are identified which allows to derive an explicit solution for 
above model. In all of these papers the dividend payments were not taken into 
account. 

We shall assume that the second company, to be called reinsurer, receives 
less premium per amount paid out, that is, 

Cl > C2. (4) 

Under the assumption made also in this paper that the premiums income 
per unit time Ci are larger than the average amount claimed A]E[J7i] the two- 
dimensional surplus in the model ([3]) has the unrealistic property that it con- 
verges to infinity on each co-ordinate with probability one. In answer to this ob- 
jection De Finetti [20J introduced the dividend barrier model for one-dimensional 
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model ((J), in which all surpluses above a given level are transferred to a ben- 
eficiary. Further, usually the payment of dividends should be made in such a 
way as to optimize the expected net present value of the total income of the 
shareholders from time zero until ruin. Associated to each dividend payment 
could be a fixed cost of size K > 0. Then it is no longer feasible to pay out 
dividends at a certain rate. Therefore the impulse control is considered instead 
where a dividend strategy is given by the pairs (Tj, Ji) for the stopping times 
Ti representing the times at which a dividend payment is made and for Ji being 
positive random variables representing the sizes of the dividend payments. 

In the mathematical finance and actuarial literature there is a good deal of 
work on dividend barrier models and the problem of finding an optimal policy 
for paying out dividends. Gerber and Shiu [23], Grandits et al. [26] and Jean- 
blanc and Shiryaev |29j consider the optimal dividend problem in a Brownian 
setting. Irback [25, Zhou [43|, Zajic [4^, Avram et al. [5], Kyprianou and 
Palmowski [31] and Loeffen [S^ study the constant barrier model for a classical 
and spectrally negative Levy risk process. Azcue and Muler [10] follow a vis- 
cosity approach to investigate optimal reinsurance and dividend policies in the 
Cramer-Lundberg model. Since under De Finetti's objective ruin is sure and its 
"severity" ignored, several alternative objectives have been proposed recently, 
involving a penalty at ruin, based on a function of the severity of ruin. This 
kind of problems is considered in Loeffen and Renaud [3S| . 

The impulse control literature is also vast. An important type of strategy 
for the impulse control problems is so-called (a+,a~) policy which is similar 
to the well known (s, S) policy appearing in inventory control models. The 
(a~, a+) policy is the strategy where each time the reserves are above a certain 
level , a dividend payment is made which brings the reserves down to another 
level a~ and where no dividends are paid out when the reserves are below a~ . 
Similar model of paying impulse dividends is considered paper where regulation 
of the risk process brings it to some fixed point. In case when risk process is a 
Brownian motion plus drift, Jeanblanc and Shiryaev [55] shows that an optimal 
strategy for the impulse control problem is formed by a (a~, a"*") policy. Paulsen 
[55] considers the case when risk process is modeled by a diffusion process and 
showed that under certain conditions a (a^,a+) policy is optimal. Note that 
in Paulsen [38] this type of strategy is referred to as a lump sum dividend 
barrier strategy. Similar model is considered in Alvarez [1], Cadenillas et al. 
[13] . Loeffen 04] models reserves by a spectrally negative Levy process and 
find conditions under which the barrier strategy is optimal. Avram et al. [9] 
find optimal impulse strategy for a Levy risk process and polynomial penalty 
function. 

Wc will follow the same arguments in the two-dimensional set-up introducing 
the following controlled risk process: 

m = {Yi{t),Y2{t)) = x{t)^m, (5) 

where IJt) = (Li(i), L2{t)) and Li{t) (i — 1, 2) are nondecreasing jFt-adaptable 
processes for a natural filtration {^t}{t>o} of 
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In this paper analogously to the one-dimensional case we will consider two 
controlling mechanism: refraction and impulse payments. In Section [2] we deal 
with the refracted process Y_ for which: 



(6) 



describes two-dimensional linear drift at rate 

^-(<^i,^2)>(0,0) (7) 

which is subtracted from the increments of the risk process whenever it enter 
the fixed set: 

B — {{x,y) : x,y > and y>b~ax}, a,b>0. (8) 

The case S_ — c — a for c = (ci, C2) and a= (—1, a) corresponds to the reflecting 
the risk process at the line y = h— ax. Then the risk process starting in U So , 
where Bq = {{x,y) '■ x,y > Q and y = b — ax}, stays there up to the ruin 
time (see Figure 1). We will also assume that the first insurance company (that 
one with lesser initial reserves) will get eventually ruined, that is: 

ci-di< 0. (9) 

Note that this approach takes into account the dependence between the reserves 
of two companies in a procedure of paying the dividends. In particular, we allow 
of moving some reserves from one company to another in the form of dividend 
payment (or moving reserves between two branches of one insurance company). 
Besides it, in the case of independent risks the barrier is usually a quadrant 
in contrast to the model considered here where a line barrier forces true join 
dividend policy of two insurance companies which are related with each other 
via proportional reinsurance. 

In the Section [3] we consider impulse controlling where Li{t) {i = 1,2) are 
cumulative payments made whenever two-dimensional risk process hits horizon- 
tal line y ~ U2. The size of the ith payment equals Ji = Xi{Ti) — ui + cie^*^ 

where Ti is the ith moment of hitting line y = M2 by the process 2£.it) e^*^ 
is independent of X exponential random variable with intensity A. Each pay- 
ment corresponds to reducing reserves to some fixed levels (-^1,1*2) and paying 
dividends while waiting for the next claim to arrive. Moreover, by condition ([4]) 
each time the second company (or first branch of company) has at least U2 of 
reserves, the first company's reserves are greater than ui. The payments of the 
dividends are always made by the first company which has greater premium rate 
(by reducing the reserves to level ui). Associated to each dividend payment is 
a fixed cost of size K > Q. In this model the reserves of the second branch of 
the insurance company serves as a random control mechanism giving us infor- 
mation when we should pay dividends from the reserves of the first branch. It 
is very convenient policy of paying dividend for insurance companies that have 
branches related via proportional reinsurance. 
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In this paper we focus on finding the 71th moment 



(10) 



of the discount cumulative dividend payments 



D= (1,1). / e 



(11) 



made until the ruin time a = ini {t > : min (li(i), l2(i)) < 0} where q is a 
discounted factor. Above (ai,Q;2) • (03,014) = 0103 + 0204 denotes the scalar 
product of two vectors in the plane. Note that a is a first time when the 
controlled risk process exit the positive quadrant and it corresponds to the first 
time when at least one company get ruined. 

The paper is organized as follows. In Section [2] we derive a partial integro- 
differential for Vn- In Section [3] we provide with explicit unique solution in 
terms of an infinite series in the case when claim sizes Ui are exponentially 
distributed. Finally, in Section|3]we solve the second dividend problem with the 
impulse control. Section [5] includes conclusions. 

2 Barrier control 

In this section for the risk process ^ we consider refracted process ([S])-® for 
the barrier set ([5]). In this section wc will derive partial integro-differential 
equation for Vn defined in (fTO|) . 

Denote u = (mi, 1*2), So — Si + S2 and 



for a general function /. 

Theorem 1 The partial derivatives ^^{u) ft = 1,2^ exist and for u S B'^ the 
function Vn is a unique solution of the equation: 





with the boundary conditions: 



nSo Vn-i{u) 




du 



ueB, 



(13) 



lim Vn{u) = 0, 



u<eB\ 



(14) 



K(0,6) = 0. 



(15) 
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Proof For y > define the Laplace transform: 

M{u, y) = E [e-y^\XiO) = u] ■ 
By the Strong Markov Property we have that for u G B'^: 
Miu, y) = (1 - \dt)M{u + cdt, ye""''**) 

+\dt / M{u + cdt~{l,l)v,ye~'''^*)dF{v) 



+Xdt / dF{v) + o{dt). 

J n\\n{u\+ c\dt ^u^+ C'^dt) 

Then the Taylor expansion and collection of terms of order dt yields 

dM dM /.min(ui,«2) 

c- ^{u.y)-\M{u,y)~qy—{u,y) + \ j M{u~ {l,l)v,y) dF{v) 

+A(1 - F(min(ui, U2))) = 0. (16) 
Moreover, for this equation we have the following boundary condition for u G S: 
M{u, y) = (1 - \dt)e-y^°'^'M{u + (c - 5)dt, ye-"''') 

nmh-i{ui-\-{ci —5i)dt,U2-{-{c2—S2)dt) 

+xdt / e-y^'"'' 

Jo 

M{u - (c - 5)dt ~ (1, l)v, ye-"''') dF{v) 



+Xdt e-y^'"'' / dF{v) + o{dt) 

J min(ui+(ci— i5i)(it,«2 + (c2— <52)dt) 



which implies: 



{c-5)-—{u,y) 



dM 

- {ySo + X)M{u,y) - qy-^{u,y) 



ueB 



+A / M{u~{l,l)v,y)dF{v) + X{l-F{mm{ui,U2)))^0. (17) 

Jo 

Note that M(u,y) and its partial derivatives with respect Ui (i — 1,2) and y 
are continuous. Thus from (|16p we finally derive for uCzB: 

(18) 

uGB 



y6oMiu,y) = 6^{u,y) 



Using representation M{u,y) = 1 + ''~n\ ^nilL) and equating the coeffi- 

cient of (— y)" in lini) completes the proof of (ll2p - (fT5|) . Moreover, conditions 
(fT4|) - p^ follows straightforward from definition of D in (fTTj) and inequality ([9]). 

We will prove the uniqueness of the solution using arguments similar like in 
Gerber [22]. It suffices to prove the uniqueness of the solution of (fT6|) and (fT8|) . 
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For that purpose, we define an operator A by 

fT />min{ui cit ,U2 + C2t) 



Jo Jo 



JO 
-At 

T 



+\E I e-^' exp <j -~ySo / e-"^' ds ^ ^ 

l(u+{c-5){t~T) ~ il,l)v,ye'~'^^)dF{v)dt, 



min(tii+(ci-5i)(t-T),ti2+(c2-<52)(t-T)) 



where T = {h — aui — u-2)j{ci + aci) for m e S'^ is a first time of getting to the 
hnear barrier 

'Bo = {(a;, y) : x, y > and y = h — ax\. (19) 

Note that in A each increment can be interpreted as a conditioning on whether a 
claim occurs before the surplus process hits the barrier (t < T) or after this event 

(in which case we have an additional term exp|— j/Jq fi""^" ds | representing 

the discounted dividends paid until the claim occurs). The solution AI of (I16p 
and (|18p is a fixed point of the integral operator A. For two functions gi and 
(?2 we have: 

1-451 (u, y) - Ag2{u,y)\ 

< llSifey) -52(u,2/)||oo I A / e^^^dt 



-At , 



+A / e-^' exp <^ ^ySo / e""" ds } dt 



T 



< \\9iiu,y) - 52 (w, 



where || • ||oo is the supremum norm over u g and y G M+. Thus it follows 
that is a contraction and the fixed point is unique by Banach's theorem. □ 

Remark 1 Equation (fT6| (hence also equation (fT2|) ) could be also derived using 
Feynman-Kac formula. 



3 An explicit solution for the exponential claims 
size and reflection 

In this section we will find unique solution of (fT2ll -(ll5 |) for n = 1, exponential 
claim size F{v) = 1 — e~"" and reflecting at the line y = b—ax, that is 5 = c — a, 
where a = (—1, a) (see Figure 1). Note that by ([7]) we have 

C2 > a. (20) 

Without loss of generality we will assume that ui < U2- In the case when 
Ml > U2 the solution of (fT2|) -(fT3 l) could be modified in the obvious way. 
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Figure 1: Controlled two-dimensional risk process. 



Let 

-[m{aci -q- X)+ac2-q- X] + ^/ (m) 

72,o(m) = o , X ; \ r ' > 0, (21) 

2((m2 + m)ci +C2+ mc2) 

where 

^72,o('^) = [m{aci—q— X)+ac2—q— X]^ +iaq[{m'^ +rn)ci+C2+mc2)] > 0. (22) 
Moreover, for 

a' = ^ < (23) 

ci + 1 ^ ' 

denote 72,0 = 72,0 (a), 72,0 = 72,o(a') and 

71,0 = 072,0 > 0, 73,fc - a72,fe = 7i,fe+i - a72,fe+i (24) 

and 

7i,o = a'72,0 > 0> 73,fe-a72,fe=7i,fc+i-'^')'2,fe+i- (25) 
We also introduce recursively coefficients: 

Do = e--.o(a)6 (£l±i)±(£L_;^^ (26) 

7i,o(ci + 1) + 72,o(a)(c2 - a) 

^ _ 73,fc+72,fc + Q 73,fc(ci + 1) + 72,fc(c2 - a) ^72..V~72..+ib /o7^ 
7i,fe +72,fe + Q! 7i,fc+i(ci + 1) -|-72,fe+i(c2 - a) 

Similarly, Dq = 1 and 

/ 73fe + 72fc + " 73 fr(ci + 1) + 72 {-(22 - a) , a v , . 
^ 7U+72,fc+" 7U+i(ci + l)+72,fc+i(c2-a) 



Theorem 2 For (ui,U2) € 'S'^ we have: 

VAux,U2) = YDk( eT'^-""^ - " e^^ '°"^^ e^^-""^ (29) 

V 7i,fe+72,fc + a ; ^ ^ 

^ , f ' 73fr+72fr + a' \ ' 

6^'''="' - e'^3..«i 6^=^.'="% (30) 

fe=o V 7i,fe+72,fe + a J 
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where 

Coo \//oo II \ 

(31) 

Proof Taking Vi :— V the equation can be rewritten as follows: 

(91/ ^min(«i,U2) 

Ci^+C2T^ {X + q)V{u) + X / y(u-(l,l)w)ae-"''d?; = 0. (32) 



du2 

Applying an operator + to the equation ([5^ gives 



OF 9y 51/ d 

L"5-2 + "=2-^ + (ci + C2)^ — ^ h (aci - g - A) — 

ouf 0U2 auiau2 aui 

d 



+ {ac2 -q- A) — l/(ii) - aqV{u) = (33) 
au2 



since 



V{ui — v,U2 — v)ae dv 

min(iii ,U2) 

V{ui — v,U2 — v)ae'~'''" dv + aV{u). 

Since Vo{u) = 1 the boundary condition translates into 

51^ 



91/ 

(ci + l) + (c2-a)= (ci + 1)t^ 



+ (c2 - a) 



du2 



(34) 



ueBo 



where is defined p9| . We will look first for a solution of equation (|33p in a 
following form: 

( Y^g7i«i _,_ Tse^^"! ) 6^=^"= (35) 

for some constants Ti and T2. Putting (|55)l into shows that 71 and 73 are 
two real roots of the equation (later we prove that they always exist): 

Ci7^ + C27I + (ci + 02)772 + (aci — q — A)7 + [ac2 — 9 — A)72 — aq = Q. (36) 

Recall that ui < U2- Then from ([5^ we obtain: 

Tl f Ci7i + C272 ~{\ + S) + — ) 6-^1"!+^-"^ 

V 71 +72 + a/ 

^^1" e-(72+a)ui+72u. ^ T2fci73 + C272 - (A + (5) 

7i + 72 + a V 



Ag73"l+72M2 __^f^^_g— (72+a)ui+72«2 _ Q 

/ + + a 



73 + 72 + a / 73 + 72 + a 
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Note that ci7i + C272 - (A + 5) + ^7^^ = C173 + C272 - (A + (5) + ^^^^^ = 0. 
Thus 

^ -Ti(73 + 72 + a) 
i 2 — 

7i + 72 + a 

and (|35p can be rewritten in the following form: 

V 71 + 72 + a / 

To fulfill boundary condition we take V = Vi which is given by ([291) -([301). It 
is not a priori clear that the solution is of the above form, but if such a solution 
can be found, we will solve ([55]) by uniqueness demonstrated in Theorem [TJ 

We will choose 71,^, 73, ^ < 71, ^ and ^[j,, 73 ^ < 7i ^ that solve the following 
equations: 

Ci7^ + C272,fe + (ci + C2)772,fe + (aci - g - A)7 + {ac2 - q - A)72,fc - aq = (38) 
and 

ci7^ + C272,fc + (ci + C2)772,fc + ("Ci - g - A)7 + (acs - g - A)72_fc - = 0, (39) 
respectively. To satisfy the boundary condition p4)) for Vi we take 



^ i?,e^^>'=(''-™i) (e^i>'="i (7i,fc(ci + 1) + 72,fe(c2 - a)) 

(73,fc(ci + l) + 72,fc(c2-a))e^^''-"^ 



fc=0 

73, fc + 72, fc + a 



and 



71, fc + 72, fc + a 

= (ci + 1) + (C2 - a) (40) 



CXD 

E D'.e''^^''^'-'^-^^ (e^U«i (7;_,(ci + 1) + 7^_,(c2 - a)) 

fc=0 

73,fc + 72,fc 



-(73,fc(ci + l) + 72,fc(c2-a))e^^''="M 

! / 



^'l.k + 72,fe + « 

= 0. (41) 

Hence ([Ml)- (1211) holds. Note that from ^ it follows that 72,^ > and 73 ^ > 
(fc = 0, 1, . . .). The coefficient E we choose in such way to satisfy last boundary 
condition ([15]). Note also that there will be no more functions being linear 
combination of functions of type (|37p such that (|40| and (|4T|) will be satisfied 
since all coefficients are uniquely determined. 

The following facts proved in the Appendix could be collected. 

Lemma 3.1 Let fc e N U {0}. Then 
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(i) 72,0 > and 7^ q > 0; 

(a) 7i,fe and 7,' ,, (i = 1,2,3, k ^ 0,1,2, ... ) satisfying and if^) . 

\39jl . respectively, always exist; 

(Hi) 72,fc+i > 72, fc > 0, 71, fc > 73, fc and 73,^+1 < 73, < 0; similarly J2,k+i > 

72,fc > 0. ilk > 73,fc o-nd y^j^^-^ < y^j^ < 0. 

To find V one needs to calculate 71,0, 72,0 a-nd Dq first. Then based on 72,0 one 
should find the smallest negative solution 73^0 of the square equation ([35]). In 
the next recursive step using (|24p by the representation 7i,fc+i = 73,/; — 072, fe + 
a72,fe+i one can identify the biggest positive 72,^+1 satisfying ([55)1 . Then ([55]) 
produces 71,^+1 and 73^fc+i. The constants Dk are given recursively by (jST]) . 
Similarly one can identify 7 ■ {{i = 1,2,3, k ^ 0, 1,2,...). In the last step we 
calculate the constant E. 

We will prove that now that series in are convergent. To do that 

we first give the following properties which will be proved in the Appendix. 



Lemma 3.2 Let k eNU {0}. Then 

(i) limfe^oo72,fc = limfc^oo72,fc = 0° onrf limfe^oo 7j,fc = limfe^oo 7^- = -0° 
for i = 1,3; 

(it) hm.^oo ^ = linifc^o, 4ff = > 1 /or z - 1, 2, 3; 

/"mi) limi.^^, = hmi,^„o = — 1 an(ilimi,^„o — \imi,^^. = — £2., 

Using the d'Alembert criterion it suffices to prove now that the series J^kLo D^e^'^-'^'^'^e^^-'''^^ 
and Y.kLo ^fee''2,fc"2 ^i^+22i±£e''3''="i given in (gSl) are convergent. Convergence 
of series appearing in pop we prove in the same way. Note that indeed, 

:^^^i+lg72,fc+lll2g-72,fcM2g7l,fc + l"lg-7l,fc"l 

Dk 

l3.k + l2.k + Oi 73,fc(ci + 1) + 72,fe(c2 - a) 



71, fc +72,fc + a 7i,fc+i(ci + 1) + 72,fe+i(c2 - a) 

g(72,fc + l-72,0(«2-&)g(7l.;= + l-7l,fc)«l ^0 as fc ^ 00, 

since by Lemma 13.21 and by observation that b> U2, 

lijn 73, fc + 72, fc + a _ g(72,fc + l-72,fc)(«2-;)) _ JJj^ g(7l.fc+l-7l,fc)"l _ Q 

fc^oo 71 fe + 72, fc + a k^co fc^oo 

and 

73.fc(ci + 1) + 72,fc(c2 - a) 
lim — ■ = const. 

fc_oo 71 fc_,_i(ci + 1) + 72,fe+l(C2 - a) 
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Similarly, 



-Dfc+i 
Dv. 



.g72,fe + lM2g-72,fe«2 g73,fc + l«lg-73,fe«l 

73,fc+i + 72,fc+i + 7i,fc + 72,fc + g 
7i,fc+i + 72,fc+i + a 73,fe + 72,fe + a 
_ 73,/c+i + 72,fc+i + a 73,fc(ci + 1) + 72,fc(c2 - a) 



(42) 



7i,fc+i + 72,fc+i + a 7i,fc+i(ci + 1) + 72,fe+i(c2 - a) 

g(72,fe + l-72,fe)(«2-6)g(73,fe+l-73,fe)Ml ^ g g^g ^ 



■ CXD. 



□ 

For the numerical analysis we assume a = 2, Ci = 4, C2 = 3, A = 1, g = 0.1. 
The values of expected dividend payments Vi(wi,U2) for ui = 1, U2 = 2 and 
Ml = 2, U2 = 3 depending on a and b are given in the Tables 1 and 2 below. 
Note that always there exists optimal choice of the linear barrier (choice of its 
upper left end (0,6) and it slope a). This choice depends on the initial reserves 
{ui,U2)- For (ui,U2) = (1,2) the optimal barrier is determined by 6 = 14 and 
a = 0.1 and for (mi, M2) = (2, 3) the optimal barrier is determined by 6 = 15 and 
a = 0.1. This is contrast to the onc-dimcnsional case where the choice of the 
barrier is given only via the premium rate and the distribution of the arriving 
claims. 



a 


b 


6 


8 


14 


15 


20 


28 


0.1 


19.85 


27.20 


34.95 


34.93 


32.48 


25.89 


0.2 


16.33 


24.31 


33.82 


34.19 


33.32 


28.03 


0.5 


11.76 


17.74 


28.98 


30.01 


32.54 


31.21 


1 


7.22 


11.40 


21.35 


22.59 


27.17 


30.07 



Table 1: Expected value of the dividend payments depending on a and b for 

fixed (^1,^2) = (1,2). 



a 


b 


6 


8 


14 


15 


20 


28 


0.1 


19.07 


27.42 


36.51 


36.58 


34.21 


27.34 


0.2 


17.17 


24.34 


35.22 


35.69 


35.01 


29.55 


0.5 


10.94 


17.50 


29.93 


31.07 


33.99 


32.78 


1 


6.59 


11.07 


21.86 


23.21 


28.19 


31.43 



Table 2: Expected value of the dividend payments depending on a and b for 

fixed (ui,U2) = (2,3). 

The Table 3 gives V{ui, U2) for fixed a = 0.9 and b = 1.8 when ui < 1*2 and 
U2 <b — aui . From this Table if follows for example that for a given slope a of 
the linear barrier it is optimal to locate around some line. In this case it is a 
line U2 = ui. 
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Figure 2: Impulse control. 



Ui 


U2 


0.2 


0.4 


0.6 


0.8 


0.9 


1.2 





2.09 


1.58 


1.11 


0.69 


0.49 


0.03 


0.1 


2.35 


1.81 


1.31 


0.86 


0.65 


0.13 


0.2 




2.06 


1.53 


1.09 


0.82 


0.25 


0.4 






1.98 


1.45 


1.20 


0.53 


0.7 








2.11 


1.83 




0.8 










2.07 





Table 3: Expected value of the dividend payments depending on ui and U2 for 

fixed a and b. 

4 Impulse control 

In this section we consider impulse controlling. The size of the ith payment 
made at epoch Ti equals Ji — Xi{Ti) — ui + cie^*'' where Ti is the ith moment 
of hitting of X_{t) in the horizontal line y = U2 and e^*'' is independent of X 
exponential random variable with intensity A. Each payment of the dividend 
corresponds to reducing reserves to some fixed levels (ui, M2) (see Figure 2) and 
paying out dividends while waiting for the next claim to arrive. By ^ at time 
Ti the first company has more than ui of reserves. Hence the payments of the 
dividends are always made by the first company which has greater premium rate 
(by reducing the reserves to the level ui). Associated to each dividend payment 
is a fixed cost of size K > Q. 

We will consider two cases, when ui > U2 and when ui < U2. In the first 
case the ruin can be only achieved by the second insurance company in contrast 
to the second case when both companies might get ruined. 

4.1 Case ui > U2 

Define: 

T+ = mf{t > : C2i - S{t) = x}, 

tZu,+. - inf{t > : C2i - S{t) < -(«2 - x)}. 
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Let a; be a size of the first claim Ui chosen according to the distribution function 
F. If there is no ruin at the moment of the arrival of the first claim then 
X < U2. When the first claim arrives the risk process jumps then from (wi,M2) 
to {ui — x,U2 — x). Note that the first portion A of the payed dividends equals 
to the discounted independent exponential random variable e^^'' = e\ with the 
parameter A multiplied by ci plus discounted jump on the line y — U2 which 
equals (ci — C2)t+ minus the costs of one payment made at the moment of 
impulse payment (that is at the moment of jump at the line y — W2). Thus, 



A = ciE 



E 



Cl 



\{c, C2)Tt - K) l(.+ <._-^^^j] dF{x) 



g + A q + \ Jo 



E 



((ci - C2)r+ - K) e-^^^ 1, 



dF{x). 



and mean of the cumulative discounted dividends equals: 

U2) ^A + A E[e-«(^-+^- ) lfr+<r- )] dF{2 



+Ap^ + Ap^ 



1=1 



A 



where 



A 



P 



q + X 

From Kyprianou 30J it follows that 



fne-'^^ \r.,.-_^^JdF{x) 

W^lHu2 - X) 



E[T+e-«^=- 1 



W(i){u2) ' 

d fW(i^U2-x) 



(43) 



(44) 



where W^?) • [q, 00) ^ [0, 00) is a scale function, that is continuous and increas- 
ing function with the Laplace transform 



> $(5), 



(45) 



for ^(6*) = 026* + \{Ee^^^ — 1) being the Laplace exponent of Xi and for 
$(g) = sup{6' > : "il^iO) = q} being its right inverse. 

For the exponentially distributed claim sizes with parameter a we have that 
tlj{9) = C2O — X9/ {a + 6) and the scale function ly^?) jg given by 



14 



where A± = 



- "+''^^^^ and 



-9- (9) 



Then, 



and 



Cl 



g + A 



C2(g + A)V^(9)M q+{q)-q-{q) 
X 



q + X 



(ci - C2) 



E 







with 



E 



dF{x) 

2 / W(«)(w2-a;)ae-"^da; 
Jo 



(W^W(U2))2C2 -<?-(<?) 



W^(9)(W2) ic2{q+{q)+a) 



C2{q 


-{q) + a) 




a 




-iq)+a) 




a 



) 



" (9)1*2 _ g-au2 



where 

q^iq)' 



1 

2^ 



C2iq~{q)+a) 

q + X + ac2 



+ {qyA+ (e«''(«)"= - e"""^) 



'(9)"2 _ g-au2 



1± 



a/((7 + A - Q!C2)^ + 4c2qa 



= fq^- = 

4.2 Case ui < U2 

Let a; be a size of the first claim U\ chosen according to the d.f. F. If there is 
no ruin at the moment of the arrival of the first claim then x <ui. Define 

Tu = mi{t > : Z{t) = U2} 
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Figure 3: Piecewise lower barrier when ui < U2. 



tl = inf{t > : Z{t) < niax(0, U2 — wi — (ci — C2)t)}, 

where Z{t) = {u2 - x) + C2t - S{t). Note that when Xi(0) = ui - x,X2(0) = 
Z{Q) — U2 ~ X > U2 ~ ui, then {Xi{t) > 0,X2{t) > 0} is equivalent to the 
requirement that {tl > t} (see Figure 3). Moreover, the impulse payment 
equals (ci — C2)tu. Thus similarly like in the previous case: 

Vi{ui,U2) = (46) 
I — p 



where 



and 



with 



p = —— / E [e"^^- dFix) 

q + Jo 



A=-^^-Kp+ ^^f^ fj E [r^e-^- dF{x) 



E [rue-'^^" l(..<r.)] = -^E [g-^^- • (47) 

To compute E [e^''^" l(T[/<ri,)] we introduce a new probability measure P'^^i^: 



dP 



^ g<i'(g)(Z(t)-Z(0))-gt^ 



On the new probability space (fi, JF, {^t}t>o, P*'"^') we have Z{t) = (m2 — x) + 
C2t-S^i\t) for 



AT,- 



and where A^**-^^ is a Poisson process with intensity Ag = A-F'(<f>(g)) for F(0) = 
/o°° e.~^^F{dx) and has a distribution function: 

Fq{dx) = e-*(«)="F(da;)/#($((j)); 

see Th. 4.8, p. 38 of Asmussen [4 and Rolski et al. [40] for details. Denote 
r,; = inf{t > : X,(t) < 0} (i = 1, 2). 
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Lemma 4.1 We have, 



where for y > U2 — ui and R — — — — 



C1-C2 ' 



P*(«)(t2 = oo|X2(0) = z)dz. (49) 
Proof Note that from Optional Stopping Theorem 

where we use fact that Z{t/\Tu) < U2. Introduce V'-'^^y) = P'^'^'^^tl = 
oo|Z(0) = y). Using strong Markov property 

since P*(«)(tc/ > tl, tl = cx)) = 0. This completes the proof of The 
identity ([iS]) is also a consequence of Markov property applied to Z at time i? 
being the zero of the line y = {u2 — Wi) — (ci — C2)t. □ 
The terms appearing in could be identified in the following way. Indeed, 
the ruin probability P*(«)(r2 < oo|X2(0) = z) = 1 - P*(9)(t2 = cx)|X2(0) = z) 
was analyzed in many papers. The reader is referred to the books by Gerber 
[21], Grandell [25j, Asmussen [4j and Rolski et al. [40]. In particular, assume 
that the distribution function F is absolutely continuous with respect to the 
Lebesgue's measure. Define gj{x) = Cj^^^^jj^ for j = 1, 2 and for x > 0, 

^* ^ d^ ^ ^ for J = 1, 2. (50) 

I— 1 

From Beekman-Pollaczeck-Khinchine formula we have 

P*(^)(r2 < (^1^2(0) = z) = fl-^ rydF,{y)) r /,^_./,,(y)e^(^/^-^)d2/; 

see e.g. Cor. 3.8 of Lefevre and Loisel [32]. For the exponential claim size with 
parameter a, Ag = Xa/{a + <i>(q)) and Fq{dx) = aqe~"''''dx for aq = a + ^{q) 
and ^q) = (A+9-^2^)+V(A+g-c2a)^+4c.,. ^ Moreover, 

P*(«)(t2 < oo|X2(0) = z) = ^ exp{-(a, - Xq/c2)z} ; 

C2ag 

see e.g. Cor. 3.2, p. 63 of Asmussen [4] and Th. 9.1, p. 108 of Asmussen [3]. 

Similarly, from the ballot theorem (see Borovkov [11], Picard and Lefevre 
[39] and Lem. 3.3 of Lefevre and Loisel [32]) we have the following lemma. 
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Lemma 4.2 For (p{z) — v ~ + R and v 



{y - {u2 



ui))/ci, 



dz 



■*(«)(ti > R,Xi{R) e dz\XiiO) = civ) 



J„ Ci(i? + W-W) 



w)fl,^y{w)dw. 



(51) 



To summarize, to find the cumulative dividend payments in this case we use 
identity ([46]) which is based on finding E [e"'''^'^ ^{tu<tl)\ identified in Lemma 
liUvia P*(«)(r2 = oo|X2(0) = z) and /'*(«)(ri > R,Xi{R) € dz\Xi{0) = c^v) 
for some v. This both quantities could be found using function // {j = 1,2) 
given in ([50]) . Unfortunately, more explicit expression for the expected value 
of the dividend payments is impossible to derive in contrast to the case when 
Ml > U2. Although De Vylder and Goovaerts [151 US] show how to analyse 
numerically expressions of type ([5T|). calculating expected dividend payments 
Vi{ui,U2) for ui < U2 still comes up against problems with identifying (|47|) and 

5 Conclusions 

In this paper we analyzed the joint dividend payments for the two-dimensional 
risk process, where the reserves of two insurance companies are related with 
each other through the proportional reinsurance. We considered two control 
mechanisms: refracting at the linear barrier and the impulse control. 

In the first case we derived partial differential equations for the nth moment 
of the cumulative dividend payments. We found also an explicit expression for 
the expected value of dividend payments when arriving claims have exponential 
distribution and risk process is refiected at the linear barrier. It appeared that 
in contrast to one-dimensional case this quantity is given by more complex 
expression being the linear combination of two series. Besides it, the numerical 
analysis showed that the optimal choice of the barrier (its upper left end (0, h) 
and its slope a) depends on the initial reserves of both companies. Thus the 
dependence of two companies might affect the choice of the optimal strategy. 

The impulse control with the positive transaction costs produces more com- 
plex expressions for the mean of the dividend payments. The expression for the 
case when the first company has less initial reserves than the second company 
appeared to be much more difficult to analyze than under the complementary 
condition. This is imposed by the way we pay dividends: the payments (reduc- 
ing reserves to the fixed level) are always made by the first company and they 
are realized when the second company has more than another fixed amount of 
reserves. Thus reserves of the second company serves us a control mechanism 
in this case. 

In future research of analyzing dividend problem for the two-dimensional 
risk process, it would be interesting to extend our analysis to the case of a more 
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sophisticated situation, such as introducing the penalty function which takes 
into account the severity of ruin, defining differently the ruin time or considering 
other types of the dependencies between reserves of the insurance companies. 
In particular, in our model the claims arrive to the both companies at the same 
time which induces jumps only in one direction. Relaxing dependence structure 
could be realized by allowing some claims to arrive to one company without 
making loss at the second company. In other words, it would be interesting 
to consider different directions of jumps. Nevertheless, we leave that point for 
the future research, since it seems that our methodology cannot be applied 
straightforward to that extended framework, and in the case of linear barrier it 
will result in more sophisticated partial differential equations. 

Appendix 

Proof o/ Lemma [XH Putting 710 = 072,0 into (|38p gives equation: 

Ii := ((a^ + a)ci + (1 + a)c2)l2.a + {aiaci - q - ^) + ctC2 - 9 - A)72,o ~ ag = 

which has one positive solution 72.0 given in (|2ip and one negative solution, 
since < 0. Moreover, by O we have 

(a^.' + a')ci + (1 + a')c2 > (52) 

and hence equation 
((a^'' + a')ci + (1 + a')c2)72,o + (^'(q^ci - q — X) + ac2 — q — A)72.o — a? = 

has also one positive solution 72 q. This completes the proof of (i). Note that 
73^0 and 73 are also well-defined since for ([38]): 

^71,0 ■= (ci - C2)^7io + 272,o(ci - C2)(aci + \ + q) + {aci - A - q)^ + Aciaq > 
and 

(ci - C2)^72;o + 272,o(ci - C2){aci + X + q) + {aci - X - qf + Aciaq > 0. 

From ([55]) we have that 

2^2 := ci7i,o73,o = 227^0 + ("^2 - 9 ~ A)72,o ~ aq <0, (53) 

where the last inequality follows from the form of Ti. Indeed, then 

2^2 = -72.o[((a^ + a)ci + 002)72,0 + a{aci - q ~ X)] 

is negative since by (i) 72,0 > and by (PT|) 

—a{aci — q — X) 
[a^ + a)ci + ac2 

= (ci — C2){aac2 + Xa^ + Xa + qa + qa^) + a{a^ + a)ci{aci — A) 
-qa{a^ + a)ci + . /A^((a^ + a)ci + 002) > 0. (54) 
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The inequality ([M]) is a consequence of the assumption the net profit con- 
dition ci > ^ and the inequality 



which is by ((22|) equivalent to 

-|-2ga(aaci + ac2 + Xa + \ + q + 2aci) + Aaqc2 > 0. 

Then by 



73,0 < 0. (55) 

By (|20|) we have a — C2 < and hence 73 q < 7i o — 0- 

Assume now that for some k E NU{0} the quantities 72, fc > 0, 73, ^ < exist. 
The solutions of equation ([55)1 then equals 

-(ci + C2)72,fc + X + q-aci- ^A^ 
73.fc = (56) 

-(ci + C2)72,fc + X + q-aci + 
7i,fc = , (57) 

where 

^■yi.k = (ci - C2)^7|fe + 272,fe(ci - C2)(aci + A + g) + (aci - A - g)^ + 4ciq:9. 

Note that 71 ^ > 73 ^ and that A^^ ^ > 0. Moreover, using (|24p we find that 

71, fc+i = 73, fc — 072, fc + a72,fc+i- Putting it into ([55]) produces the equation for 

72, fc+i: 

((a2+a)ci + (l + a)c2)72%+i 

+ [(73,fc - a72,fc)(2cia + ci + C2) - (A + q){l + a) + a{aci + C2)]72,fc+i 
+ (ci(73,fc - 072, fc)^ + {cia - A - g)(73,fc - 072, fc) - aq) ^ (58) 

which has a solution since: 

^72,fc+i = ["(aci + C2) - (A + g)(l + a)]2 

+4Q;g((a^ + a)ci + (1 + a)c2) + (73, fc - a72,fe)^(ci - 

+2(73,fc - a72,fc)(c2 - ci)((A + g)(l + a) + a{c2 + aci)) > 0. (59) 

Furthermore, 



72,fc+i - 72,fc = \/ ^-f2.k+i + (2cia + ci + C2)J A^,_^ 



+ ^^^if^72.fc + ^^^(a + A + ,) > 0. (60) 
2ci 2 
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Hence 72,fc+i > and 

^7i,fc+i = (ci-C2)^72,A;+i+272,fc+i(ci-C2)(aci+A+(7) + (aci-A-9)2+4cia(7 > 

which means that there exist two solutions 7i,A;+i, 73,fc+i < 7i,fc+i of equation 
([55]) . Similarly, by ([5^ and ([SO)) , using the same arguments we can prove that 
72 ^_|_]^ > 72 fc ^^"^ that that there exist two solutions j[ f.^■^^,J'^ < j[ f.^^ of 
equation ([55)) . This completes the proof of (ii). To prove (iii) note that: 



73,fc+l - 73,fc = (Cl + C2)(72,fe - 72.fe+l) + -^/a^- ^. 



A^i,,+i < 



since y^A^^ ^ ~ y^A^JTT+T < which is equivalent to the inequality 

(Ci - C2)^(72,fc - 72,fc+l)(72,fc + 72,fe+l) 

+2(72, fc - 72,fe+i)(ci - C2)(aci + X + q) <0 

that foUows from ()60p . In the same way we prove inequality 73 fc+i < 73 fc- D 
Proof of Lemma ] 3. 2\ Recall that by Lemma IS.lT iii) the sequence {72. fc,^ = 
0, 1, 2, . . . } is monotone. Hence if it there were not true that limfc^oo 72_fc = 00, 
then it would exist C > such that C = limfe_»oo 72, fc- Then by ([55)) . 



lim (72,fc+i -72,fc) = lim ( JA^^.^+i + (2cia + ci + C2) JA^.^ J 

, (Cl - C2)^ (C1-C2) 

H — C H [a + A + g) > 0. 

But if hmfc^oo72,fc = const then limfc^oo(72,fc+i - 72. fc) = which gives a 
contradiction. Note that by ([5S)) we have that limfc^oo = ~1 and hence 
limfc^oo 73,fc = -00 and by ([58|), 



lim 1M±1 = ^£l±^ > 1. (61) 

fc— 00 72^fe aci + C2 

Similarly, by ([57]) limfe-00 ^ ^ -ff, hence limfc^oo7i,fc = -oo- Moreover, 
hmfc^oo ^^^i^ = -^2£i±£i and thus hmfe-00 ^^^^^ = cia±£i ^ Similarly we 

72, fc Cl aci+C2 J. cia+C2 

can prove that limt^oo '^^''''^^ — £l2±£i_ xhe same limits we can obtain by 

jr- — JO Cia + C2 

exchanging ji^k by 7^' j, for i = 1, 2, 3. □ 
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